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We investigate whether the quantal and experimental am- 
plitudes of current oscillations of Resonant Tunneling Diodes 
in tilted fields are obtainable from Periodic Orbit (PO) the- 
ories by considering recently proposed PO approaches. We 
show, for the first time, that accurate amplitudes and fre- 
quency shifts for the current oscillations (typically to within 
a few %) can be obtained from a simple analytical formula 
both in the stable (torus-quantization) limit and the unstable 
regimes of the experiments which are dominated by isolated 
PO's. But we find that the PO approach does not describe 
quantitatively the dynamically interesting intermediate ex- 
perimental regimes which appear to be dominated by contri- 
butions from complex orbits and multiple non- isolated PO's. 
We conclude that these regimes will not easily be described 
by the usual PO approach, even with simple normal forms. 
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Mesoscopic systems have been extensively investigated 
as probes of 'quantum chaos' in real systems. Among 
these, the Resonant Tunneling Diode (RTD) in tilted 
fields, first introduced in 1994 H , has attracted much 
attention because of the diversity of observed effects it 
exhibits which have been attributed to Periodic Orbits 
and 'soft chaos'. For instance, it has been used as an 
experimental probe of spectral fluctuations due to unsta- 
ble periodic orbits 0], quantum scarring ^-Q], bifur- 
cations , 'ghosts' |)],[l0) and the torus quantization 
regime |Kj ■ All these effects manifest themselves through 
observed oscillations in the tunneling current of varying 
amplitude and frequency. 

The well-known Gutzwiller Trace Formula (GTF) Q 
is a powerful tool in the quantization of chaotic systems, 
ft relates the frequencies and amplitudes of oscillations in 
the density of states to the actions and stabilities of clas- 
sical periodic orbits (PO's) in a simple analytic formula. 
However to date we have no equivalent formula which 
describes the corresponding oscillations in the tunneling 
current. Hence some interpretations of the experiments, 
for instance whether bifurcating orbits are seen or not, 
remain controversial 
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In JT3] a semiclassical 
treatment of the current was presented, evaluating the 
Bardeen tunneling matrix element numerically within the 
Wigner phase-space representation. Thus they could re- 
produce some interesting experimental features such as 
regimes of period-doubling of the current. 

We show here that a simple analytical formula devel- 



oped by Bogomolny and Rouben jl5| can yield quan- 
titatively the amplitudes of contributions from isolated 
periodic orbits - stable or unstable. This represents the 
first demonstration that the amplitudes of the current 
oscillations may be quantitatively described by a peri- 
odic orbit expansion. The formula is in good agreement 
with results obtained using the approach of |14|. The 
latter approach provides more flexibility allowing, for in- 
stance, excited initial states (which are not required in 
our calculation). But the advantage of the new formula 
lies in its simplicity -it is as easy to evaluate as the GTF- 
and in the physical insight which the analytical expres- 
sion provides. For instance, it exposes a shift between 
the frequency of the Gutzwiller density of states oscil- 
lations and the current oscillations. This shift is small 
(< 1%) but the resolution of the quantum scaled calcu- 
lations easily exceeds this. The comparison carried out 
here confirms that this shift is indeed easily detectable 
and accurately predicted in terms of classical quantities. 
Also, the formula reproduces observed features due to 
torus-quantisation, in agreement with quantal calcula- 
tions and the simple model we proposed previously [flO| . 

We have applied here a rigorous test to this formula, 
since we have compared it with a broad range of accu- 
rate amplitudes obtained from a scaled quantum spec- 
trum |9T]i"(i(|. We also compare with the extensive set of 
experimental data obtained at Bell Labs 0. 

Many of the most interesting experimental features, 
such as period-doublings occur in an intermediate regime 
characterized by contributions from multiple non-isolated 
orbits and complex PO's (ghosts). We find here that 
agreement in this regime is qualitative: both the formula 
and the approach of jl4j yield the rough range of period- 
doubling regions, but the amplitudes are in poor agree- 
ment with quantal results and experiment. In particular, 
in the two 'ghost' regions we identify, we cannot account 
for the amplitude of the current oscillations, even with 
normal form corrections. The strength and persistence 
of these contributions remains one of the most puzzling 
features of these experiments since in general 'ghosts' 
are strongly damped away from the tangent bifurcation 
where they appear. 

We recall briefly the RTD model |Q. In essence, the 
physical picture is as follows: an electric field F (along x) 
and a magnetic field B in the x~z plane (at tilt angle 6 to 
the x axis) are applied to a double barrier quantum well. 
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Electrons in a 2DEG accumulate at the first barrier and 
tunnel through both barriers giving rise to a tunneling 
current /. In the process they probe the classical trajec- 
tories -regular or chaotic- arising from specular reflection 
at the barrier walls. The current oscillates as a function 
of applied voltage V. After rescaling with respect to B, 
the dynamics at given 9 and ratio of injection energy to 
voltage (R = E/V ~ 0.15 for the Bell Lab experiments) 
depends only on the parameter e = V/LB 2 . Regular 
behaviour occurs at high e ~ 20000 (in atomic units), 
chaotic behaviour for low e ~ 1000. The well width is 
L = 1200A As we study here small angles (9 < 27°), we 
neglect the shift Sz ~ d tan 8 due to the mean distance d 
between the 2DEG and the left inner barrier (x = 0). At 
9 = 0° the current consists of pure period one oscillations, 
of amplitude independent of e, associated with a straight 
line PO (to), bouncing alternately between walls. This is 
our reference current amplitude In and we normalise all 
amplitudes (semiclassical, quantal and experimental) to 
Io. 

As 9 increases to is no longer a straight line but contin- 
ues to dominate the period-one oscillations. Due to the 
relatively short coherence time r ~ O.lps 0] just four 
of the shortest PO's (to, its second traversal 2to and the 
period-two PO's Si and S') account for all experimental 
features studied here. Their shape is shown in Fig. [l|. 
Nevertheless, their dynamical behaviour is far from sim- 
ple. Already in |lj it was observed that these classical 
PO's appear and disappear abruptly, to undergoes a se- 
ries of tangent bifurcations, where it ceases to exist as a 
real PO but leaves a complex 'ghost'. Subsequently, be- 
low the tangent bifurcation, a new similar looking PO t' 
re-appears from the opposite side of the Surface of Sec- 
tion (SOS) and re-stabilizes . S' appears abruptly at the 
discontinuity in the potential between a barrier and the 
energy surface due to the magnetic confinement (a 'cusp' 
bifurcation J^]). It disappears subsequently in a tangent 
bifurcation also leaving a 'ghost'. 

Despite these intricacies, we can identify in the ex- 
periments three generic dynamical regimes. These are 
illustrated in the Surfaces of Section in Fig. |l|. They 
are: (1) e > 20000. The torus regime: the large sta- 
ble island of to yields a period-one current. The current 
shows 'jumps' associated with torus quantization. (2) 
20000 > e > 3000. Intermediate regime, with contribu- 
tions from multiple non-isolated PO's or complex orbits. 
(3) e - 3000 - 1000. Unstable regime. We can identify 
contributions from unstable, isolated PO's such as Si 
(the period-two oscillation identified in the original Not- 
tingham experiments M). We show below that the new 
semiclassical formula gives excellent results in (1) and (3) 
but rather poor results in the intermediate regime. 

The theoretical scaled current (neglecting experimen- 
tal broadening due to incoherent processes) is a den- 
sity of states weighted by a tunneling matrix element: 
= Si Wi6(N - Ni). 



(D 



(2) 



(3) 




FIG. 1. Shape in x — z plane of the main PO's to (pe- 
riod-one), S' and 5*i (period-2). The SOS's illustrate, 
for 9 = 11°, three generic dynamical regimes typical of 
6 = 10-30°. (1) e = 20000. Large stable island. (2) 
e = 7000. Intermediate regime, to about to undergo a tangent 
bifurcation near the edge of the SOS, which will remove the 
real PO. This regime is characterised by contributions from 
non-isolated POs or complex orbits. (3) e = 3000. Unstable 
regime. The new t' PO has re-appeared on the far side of the 
SOS. A strong period-two signal from the isolated unstable 
PO Si, which occupies the central region most accessible to 
the tunneling electrons, is seen in the experiments. 

The experimental range (V = 0.1 — 1.1 V) corresponds 
to N ~ 12 — 43, which gives an average N ~ SU 1 ~ 28 
corresponding to V = 0.5 V. In fact N is a re-scaled 
magnetic field §,0: N = BLy / 2mLe(R + 1/2) /it. In 
our calculations we used the Bardeen matrix element filf | 
form for the tunneling probability. Then as explained in 
]l5| , one can re-express the matrix element in terms of en- 
ergy Green's functions and use their semiclassical expan- 
sion over classical paths. We consider the initial state de- 
scribing the electrons prior of tunneling to be the lowest 
Landau state: <fia(z) = B cos 9 / tt exp(— B cos 9z 2 /2). 
Then the tunneling current is given by: 



1(B) oc 5fte J dz! dz' 



-1/2 AS(z, z') -B cos 8(z 2 +z' 2 )/2 



c\(z — >z' ) 



where m u = ^. 

The integrals were evaluated analytically by station- 
ary phase with the condition dS/dz = dS/dz' = 0. This 
condition implies that only PO's starting with null mo- 
mentum p z = contribute. The resulting contribution to 
the normalized current I n orm(B) = 1(B) /Io for a given 
periodic orbit is approximated by: 

e iB(S+AS)+i^ir /2-B£ 

Inorm(R) — , 

\J — cos 9m\2 + m2i / cos 9 + 2im,\\ 
AS = cos 6^7/(1 + 7 2 ) ; £ = cos6z^ 2 /(l + 7 2 ) 
_ mil - 1 
cos6'mi2 ' 

where /i is a Maslov index, S and rriij are the scaled 



2 



action and element of the classical monodromy matrix 
of the PO, with starting position (x = 0, z — z ). The 
semiclassical theory predicts that the frequency of the 
current oscillations is shifted relative to the scaled action 
S by AS. 

We show in Fig. the scaled action of to, the semi- 
classical frequency S + AS and the quantum frequency 
obtained by Fourier transform. The shift of frequency is 
~ 1% at most, but clearly we can see that the shifted 
frequency is in excellent agreement with the quantal re- 
sults. 
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FIG. 2. Frequency of period one oscillation at (a) 6 = 11°, 
(b)20° and (c) 27°. Squares = quantum frequency (QM). 
Dashed line = action of the to (PO). Solid line = semiclassical 
frequency (SC), which is shifted relative to the action of the 
classical PO but is in excellent agreement with the quantum 
results. 



In 1 10 1 we explained how one may extract experimen- 
tal PO amplitudes in the case where the current has just 
a pure period one or two oscillation, by removing the 
smooth non-oscillatory component. This is only possible 
in a restricted range of the experiments. Also, to compare 
with theory we must in addition consider two factors: i) 
the experimental features are displaced to a lower voltage 
relative to the theoretical value V — FL; ii) incoherent 
processes damp each PO contribution by a factor e _T / r 
where T is the period of the PO. For the voltage dis- 
placement, we found that all PO features appearing at 
V = 0.5 V in the calculated spectra are systematically 
displaced to a voltage 30% smaller in the experiment. 
For example in || , where characteristic line profiles were 
correlated with different dynamical regimes, we showed 
that the distinctive spectral signature of a bifurcation, 
seen in the quantal spectra at e = 13000, is seen in the 
experiment at e = 10000. Much of this voltage shift is 
accounted for by the voltage dependence of the effective 
mass [[To). Hence we took e — » 1.3 x e for the experiments, 
for all angles and all values of e. Then we found that, for 
all angles 9 = 11 — 27°, the position of the theoretical 
and experimental period-doubling maxima are in good 
agreement. 



For the damping we find that the amplitudes of the 
maximum period-doubling current would coincide if we 
chose r in the range 0.10 — 0.12 ps for a given angle. This 
is remarkably consistent with the expected r ~ 0.1 ps 
suggested in (!]]. We chose a representative value r = 0.11 
ps for all the experimental amplitudes. In effect period- 
two amplitudes are damped by incoherent processes by 
about an order of magnitude relative to Iq . 
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FIG. 3. Quantal (QM), semiclassical (SC) and experimen- 
tal amplitudes for period one oscillations at 6 = 11° for (a) 
regimes (2) and (3) with tangent bifurcation (T), complex 
orbit (CO); (b) torus regime (1). The semiclassical formula 
gives excellent results where the PO is isolated ie for e > 8000 
(to) and e < 2000 (t' ) but gives poor results near the bifurca- 
tion and for e = 6500 — 3000 in the CO region, even including 
the ghost PO. Reliable experimental period one amplitudes 
are only obtainable at 0.5 V from e = 4000 up to the start of 
a period-doubling region at e ~ 9000. 

We show in Fig. || the amplitudes of the experimental, 
quantal and semiclassical period one current at 9 = 11° 
from the torus regime through to the unstable regime. 
For high e [regime (1)] the semiclassical current in Fig. 
|3|(b) is due to the stable to and its tori, and the agreement 
with the quantal calculation is excellent. However, the 
most interesting dynamical region is at e = 3000 — 6000, 
where there is a significant quantal/experimental current 
but no real to PO (to disappears in a tangent bifurcation 
at e = 6500, and reappears at e = 4300; between 4300 — 
3000 the new to is not easily accessible to the tunneling 
electrons so its contribution is negligible). Here [regime 
(2)] the semiclassical result is poor even when we include 
the 'ghost' complex PO. Even with a cubic normal form, 
which we do not present here, agreement remains poor. 
Quantitative agreement is once again good when the 're- 
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born' isolated real t' orbit dominates the current for e < 
2500 [regime(3)]. 
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FIG. 4. Quantal (QM), semiclassical (SC) and experimen- 
tal amplitudes for period two current at (a) 9 — 11°, (b) 
9 = 20° and (c) 9 = 27°. As previously we have 3 regimes. 
(1) stable and (3) unstable isolated PO, where the semiclas- 
sical formula is good. In regime (2) in contrast we have 
the non-isolated contributions of 2to and S' which cannot 
be added straightforwardly. Hence we simply show the in- 
dividual contributions. Here agreement between the formula 
and quantum calculations/ experiment is poor. S' appears 
abruptly at a cusp bifurcation (e ~ 18000) and disappears 
in a tangent bifurcation at lower e, below which the experi- 
ment shows a slowly decaying 'plateau' due to a complex orbit 
(CO). 

In Fig. H we show a comparison between experimental, 
quantal and semiclassical amplitudes of the period two 
current at 6 = 11°, 20° and 27°. We were unable to read 
reliable experimental period- two amplitudes for 6 = 11° 
as a strong period one beat is also present. At = 20° 
and 27°, the quantal calculations and the experiments are 
in very good agreement. As expected in the large stable 
island regime e > 25000 [regime (1)] agreement between 
the semiclassics and quantum is excellent. This is also 
the case for e < 3000 [regime (3)]. Here, the isolated 
unstable PO S\ describes the current very well. 

However, in the intermediate regime (2), the quan- 
tum current requires a coherent superposition of the non- 
isolated PO's 2io and S'. A straightforward sum (al- 
lowing for their phases) yields poor results. 2i and S' 



have near identical actions, unresolvable in the quantum 
Fourier transform spectroscopy. The 2to contribution 
comes from Miller tori localized on the large stable is- 
land and, due to the moderate values of h, substantially 
beyond the island boundary. A phase-space analysis with 
Wigner and Husimi functions shows that the S' scars are 
mixed in with the outer tori of the to island. Hence both 
their action and phase-space localization coincide. At 
11° however, the contribution of S' is small and occu- 
pies a narrow range in e. In this case the semiclassical 
amplitudes are quite good. This is not the case at 27°. 
Both the individual island (2to) and S' contributions are 
significant between e = 20000 — 12000 and there is no 
agreement with the quantal results. For e < 8000 the 
2f contribution is negligible and S' has disappeared into 
the complex plane in a tangent bifurcation. Even includ- 
ing the S"-ghost, we were unable to obtain quantitative 
agreement in this complex orbit (CO) region spanning 
e = 8000 - 5000. We note that the slowly declining 
'plateau' seen quantally and in the experiment is a sur- 
prising and unexpected feature, since 'ghosts' should be 
exponentially suppressed as the imaginary component of 
the action grows. 

We have also compared with the alternative semiclassi- 
cal approach in [ p"!] . All regimes considered here and in 
fuj involved PO's with initial p z = 0. In that case, and 
assuming the lowest Landau state for the initial state, we 
have found that the numerical integral (5) of |l4| reduces 
to the analytical formula here. Both encounter the same 
difficulties in regime (2). We note that although jL5j re- 
quires the p z — selection rule, Jl4j requires only p z to 
be small. As yet we have no unambiguous experimental 
detection of a PO with p z > 0. 

Finally, we note that in general, the approach of pl| 
would predict complex stationary phase points. Their 
complex part has been neglected in order that the theory 
may obtain PO's. Our work jlTj indicates that the con- 
sistent failure of the PO formalism in the intermediate 
regime, despite the usual normal form corrections, may 
require that the usual PO picture be partly abandoned 
since complex non-periodic contributions -as opposed to 
'ghosts' which are complex PO's- may be essential. 
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